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ABSTRACT 

We present preliminary results of the first and on-going radial velocity survey for circumbinary 
planets. With a novel radial velocity technique employing an iodine absorption cell we achieve an 
unprecedented RV precision of up to 2 ms^-'^ for double-lined binary stars. The high resolution spectra 
collected with the Keck I/Hires, TNG/Sarg and Shane/CAT/Hamspec telescopes/spectrographs over 
the years 2003-2008 allow us to derive RVs and compute planet detection limits for ten double-lined 
binary stars. For this initial sample of targets, we can rule out planets on dynamically stable orbits 
with masses as small as ^--^0.3 to 3 M-jup for the orbital periods of up to ~5.3 years. Even though the 
presented sample of stars is too small to make any strong conclusions, it is clear that the search for 
circumbinary planets is now technique-wise possible and eventually will provide new constraints for 
the planet formation theories. 

Subject headings: binaries: spectroscopic — stars: individual (HD9939, HD13974, HD47415, HD78418, 
HD86146, HD195987, HD210027, HD214686, HD221950, HD282975; techniques: 
radial velocities; stars: planetary systems) 



1. INTRODUCTION 

Searches for planets in close binary systems ex- 
plore the degree to which stell ar multiplicity in- 
hibits or promotes planet formation (|Muterspaughll200"5l : 
iMuterspaugh et al.ll2006l . l2007f) . Detection of giant plan- 



ets orbiting both components of short period (P < 60 
days) binaries ("circumbinary planets") will have sig- 
nificant consequences for theoretical understandings of 
how giant planets are formed. The binarity of the cen- 
tral body creates an environment in which the evolution 
of a protoplanetary disk is subs tantially d ifferen t than 
around single stars (|Artvmowicz fc Lubowl [r994[ ) . This 
must have an effect on the migration of giant planets in a 
disk as well as on the "parking" mechanism and their fi- 
nal orbit. Likely, also the dynamical interaction between 
protoplanets and then planets in a multi-planet system 
should be affected by the central body binarity and pre- 
sumably result in a different distribution of the orbital 
elements of planets. Finally, if one assumes that plane- 
tary orbits are coplanar with the orbit of an eclipsing bi- 
nary, then there is an enhanced probability of detecting a 
circumbinary transiting planet (,Konacki 2009: Ofir 2 009, 
[20M [Sch^^eldal[l99l ISchneider fc Chevretonl[l990h . 

The recent work by Konacki (2005,2009) demonstrates 
a method for obtaining radial velocity (RV) precisions 
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of up to 5ms~^ for double- lined spectroscopic binaries 
(SB2s) (now improved to 2 ms~^) and for the first time 
opens the opportunity to search for circumbinary plan- 
ets via RVs. This method has been applied to a new 
search for circumbinary planets: The Attempt To Ob- 
serve Outer-planets In Non-single-stellar Environments 
(TATOOINE). 

Planets in binary and multiple stars have been a sub- 
ject of theoretical works long befo re any extrasolar plan- 
ets were detected. iDvorakI (|1982t ) investigated dynami- 
cal stability of a wide range of planetary configurations 
in the framework of the elliptic restricted three-body 
problem. These included the so-called P-type (Planet- 
type, circumbinary orbits), S-type (Satellite-type, cir- 
cumprimary or circumsecondary orbits) and L-type or- 
bits (Librator-type, orbits around stable Lagrangian 
points L4 or L5 for the mas s ratios n < 0.04) . This 
subject was continued by e.g. iHolman fc Wiegertl ()1999l ) 
who based on extensive numerical simulations provided 
useful formulas allowing one to compute stable regions 
for among others circumbinary planets. 

Circumbinary disks and planet formation in such 
disks have been of some interest to theorists as well. 
The evolution of a c i rcumb inary disk is studied by 
lArtvmowicz fc Lubowl (119941 ) who demonstrates that 
such a disk will be truncated at its inner edge by tidal 
torques to within 1.8-2.6 times the sem imaj or axis of the 
binary . iMoriwaki fc Nak agawa* l|2004f ) and iSchoU et all 
l|2007h show that planetessimal acc retion should be pos- 
sible in circumbinary disks and Ouint ana fc Lissaued 
(|2006l ) provide numerical proofs that planetary systems 
similar to those around single stars may be formed 
around binary stars with the apastr on distance < 0.2 
AU. Most rece ntly, in a series of papers iPierens fc NelsonI 
(|2007l I2008 a'l^ investigate formation, migration and 
evolution of planets in circumbinary disks. In particu- 
lar, they suggest that circumbinary planets may be more 



2 



M. Konacki, M. Muterspaugh, S. Kulkarni and K. Helminiak 



common in the Saturn-mass regime as higher mass plan- 
ets are more Uk ely to undergo close enco unters with the 
secondary star (jPierens fc Nelsonl[2008Hl . 

There are known a few cases of young spectroscopic 
binaries with circumbinary disks. These include AK Sc o 
(lAndersc n ct al.. 1989), GW Q ri (.Mathieu et aUfTOOl 
DQ T au (jMathieu et aIIIT997l ). GG Tau (jPutrev et all 
Il994f ). The observations confirm the truncation of the 
disks at their i n ner ed ge as predicted by theory. Recently, 
IKastner et aP pOOl have detected CO, HON, CN and 
HCO+ around a 12 Myr-old close binary V4046 Sgr 
demonstrating that it is surrounded by a rich molecular 
disk and showing a great similarity to the ~ 8 Myr-old 
star TW Hya; a star/disk system regar ded a representa- 
tive o f the early solar nebula. Finally, 'Z uckerman et al] 
(|2008f l provide evidence that an SB2 BD-h20 307 is an 
old binary star surrounded by a debris disk likely formed 
in a collision involving a terrestial planet. 

Clearly, there is enough evidence that circumbinary 
planets should form, evolve and survive on stable or- 
bits around close binary stars. Yet no radial velocity 
survey has been carried out to detect such planets de- 
spite the fact that the RV technique for single stars has 
allowed for a thriving scientific endeavor over the last 
15 years. The fundamental problem with double-lined 
spectroscopic binary stars is that their spectra are highly 
variable due to the orbital motion of their components 
resulting in Doppler shifts typically up to ~100 kms~^ 
for each component spectrum. Hence, the approach used 
for single stars where a Doppler shift of an otherwise con- 
stant shape (spectrum) is found cannot be applied. 

In this paper we present the preliminary results of the 
first radial velocity survey for circumbinary planets. In 
§2 we discuss the novel iodine cell based approach that 
allows us to measure precisely RVs of SB2s. In §3 we 
describe in more detail our on-going RV effort to detect 
circumbinary planets. In §4 we show the planet detection 
limits for 10 SB2s from our survey and conclude in §5. 

2. PRECISION RVS OF DOUBLE-LINED SPECTROSCOPIC 
BINARIES 

In the iodine cell (I2) technique, the Doppler shift of a 
star spect rum AA is determined by solving the following 
equation (iMarcv fc Butledll992D . 

lobsiX) = [/.(A + AA,)T/,(A + AA/J] ® PSF (1) 

where AA^ is the shift of the star spectrum, AA/j is the 
shift of the iodine transmission function Tj^ , represents 
a convolution and PSF a spectrograph's point spread 
function. The parameters AAs,AA72 as well as param- 
eters describing the PSF are determined by performing 
a least-squares fit to the observed (through the iodine 
cell) spectrum lots- For this purpose, one also needs (1) 
a high signal to noise ratio (SNR) star spectrum taken 
without the cell Is which serves as a template for all the 
spectra observed through the cell and (2) the I2 transmis- 
sion function T/^ obtained, for example, with the Fourier 
Transform Spectrometer at the Kitt Peak National Ob- 
servatory. The Doppler shift of a star spectrum is then 
given by A A = AXg — AXj^ . Such an iodine technique can 
only be applied to single stars. This is dictated by the 
need to supply an observed template spectrum of each 
component of a target binary star in Eq. [TJ 
In the case when a composite spectrum of a binary star 



is observed, the classical approach with the iodine cell 
cannot be used since it is not possible to observationally 
obtain two separate template spectra of the binary com- 
ponents. This problem can be resolved as follows. First, 
two sequential exposures of each (binary) target are al- 
ways taken — one with and the other without the cell. 
This way one obtains an instantaneous template that is 
used to model only the adjacent exposure taken with the 
cell. The derived Doppler shift, AA^ (where i denotes the 
epoch of the observation), carries no meaning since each 
time a different template is used. Moreover, it describes a 
Doppler "shift" of a composed spectrum that is typically 
different at each epoch. However, the parameters — in 
particular the wavelength solution and the parameters 
describing PSF — are accurately determined and can be 
used to extract the star spectrum, ^*j!j(A), for each epoch 

&) = [lobsW PSF^Tj^iX), (2) 

where (g)"^ denotes deconvolution, and PSF'- represents 
the set of parameters describing PSF at the epoch i. Such 
a star spectrum has an accurate wavelength solution, is 
free of the I2 lines and the influence of a varying PSF. 

In the next step, the velocities of both components 
of a binary target are measured with the well know 
two-dimensional cross- correlation technique TODCOR 
(jZucker fc Mazehl[T99l using as templates the synthetic 
spectra derived with the ATLAS 9 and ATLAS 12 pro- 
grams (jKurucj ll995l and matched to the observed spec- 
trum, /s(A). This approach results in an RV precision of 
20-30 ms-i (|Konackill20"05l . Now, having the RVs ob- 
tained with TODCOR, one can carry out a tomographic 
disentangling of the composite spectra of a binary to 
produce the real (observed) component spectra and use 
them in the original equation [TJ This finally allows us 
to directly measure the RVs without relying on synthetic 
templates. Such a procedure results in RVs having a pre- 
cision of up to several m s~^ (jKonac kil [2009l ). The formal 
errors of the velocities are derived from the scatter be- 
tween the velocities from different echelle orders. 

3. THE TATOOINE SURVEY 

The TATOOINE survey was initiated in mi d 2003 with 
the 1 0-m Keck I and its Hires spectrograph (jVogt et al.l 
[T99I as an addition to an RV survey for planets or- 
biting components of speckle binary and multiple stars 
(jKonacki 2005b). The original sample of 12 SB2s was 
selected primarily to provide precision RVs for the as- 
trometric data collected at the Palomar Testbed Inter- 
ferometer (jColavita et al.l 119991 : iKonacki fc Land l2004f) 
and to test the new RV technique on a sample of close 
binaries. The Keck I/Hires survey was continued un- 
til mid 2007. The sample was also monitored from mid 
2006 until mid 2007 wi th the 3.6-m TNG t elescope and 
its Sarg spectrograph (|Gratton et al.|[200ll ). Since fall 
2006, the TATOOINE project has been collecting ve- 
locity measurements at Lick Observatory using the 3-m 
Shane and 0.6-m Coudc Auxiliary Telescope (CAT) tele- 
scopes and the Hamilton Spectrograph ( VQgt,1987> Ham- 
spec;). Overall, ~ 50 northern binaries are currently be- 
ing monitored by the TATOOINE project. Recently, we 
have also initiated a southern effort at the 3.9-m AAT 
telescope with the UCLES spectrograph (jDiego et al.l 
,1990) . 
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The spectrographs employed in the survey, Hires, Sarg 
and Hamspec, are all high resolution echelle spectro- 
graphs equipped with iodine absorption cells. They pro- 
vide spectra with a resolution of 67 000, 86 000 and 
60,000 respectively across a ^400-900 nm bandwidth cen- 
tered on the 500-600 nm range corresponding to an iodine 
cell's spectral features. The cells are used to superimpose 
a reference absorption spectrum in a manner described 
in §2. For each measurement, back-to-back spectra of the 
target binary were acquired, first with an iodine absorp- 
tion cell in the path of the starlight, then with the iodine 
cell removed. Depending on the telescope, different ex- 
posure times were used to balance the resulting SNR per 
collapsed pixel and the efficient use of telescope's time. 
The highest SNR of ~250 was achieved with the Keck 
I/Hires also to test the data pipeline in the regime of 
high SNR spectra. The typical SNR was --75-150 for the 
TNG/Sarg and -50-150 for the Shane/CAT/Hamspec. 
In consequence, the RVs have different precision ranging 
from ~2 (HD195987, the best case) to 20 ms'^ for the 
primary stars. In any case, the precision is sufficient to 
detect planets with masses as small as 0.3 Mjup- Note 
also that due to a brightness ratio between the primary 
and secondary, the RVs of the secondary are typically 
of lower precision as the SNRs are for the composite ob- 
served spectra. For example, an SNR of 250 and a bright- 
ness ratio of 6.7 (HD 195987) corresponds to an SNR of 
220 for the primary and only 30 for the secondary. 

A circumbinary planet will exhibit two indirect effects 
on the RV of the stars. First, the apparent system veloc- 
ity will vary in a periodic manner due to the motion of 
the binary about the syst em barycenter, with amplitude 
(|Muterspaugh et al.|[2007l ) 



Avb — 57 ms 



{Mp/Mjup)siiiip 



y/mt + Mp) /Mq) (ap/lAU) 



(3) 



Differential reflex motion and perturbations of the binary 
orbit by the planetary companion are expected to be neg- 
ligible on reasonable timescales. Second, the finite speed 
of light will cause apparent changes in the phase of the 
binary orbit due to the reflex motion of the binary about 
the system center of mass. This pha se shift is detected 
for p lanets with masses as small as (jMuterspaugh et all 
[2001 

(a„,/20m/s) (Pf,/5d)^/' (M,/Mq)'/' 



Mr, 



70Mjup X 



\/N — 6 sin ii, sin ip (op / 1 AU) 



(4) 

assuming 1% FAP (SNR=5.8) where TV is the number 
of observations and if, and ip are the inclinations of the 
binary and planet orbits, respectively. Note that this 
observable is actually more sensitive to longer period 
planets. In combination, these effects allow us to detect 
Jupiter mass planets or smaller for an extended range of 
orbital periods. In the current analysis, the second effect 
has been ignored as it is small on the relatively short 
timescales being considered. 

4. PLANET DETECTION LIMITS FOR TEN SB2S 

For this initial sample of stars, we have selected these 
SB2s from our survey that have long time span RV data 
sets and an orbital phase coverage allowing for a reli- 
able tomographic disentangling. These are HD 9939, HD 



13974, HD 47415, HD 78418, HD 86146, HD 195987, HD 
210027, HD 214686, HD 221950 and HD 282975. Let us 
note that we recently made improvements to our data 
pipeline and the Keck I/Hires RVs used in this paper are 
ty pically several ms~^ more accurate than those shown 
in lKonackH ()2009[ ). The internal RV errors are computed 
from the scatter between the echelle orders used in the re- 
duction. They are expected to underestimate the real RV 
scatter. In addition to the stellar jitter, one of the main 
reasons for the underestimation of the errors are the im- 
perfections in the disentangled template spectra. While 
these imperfections are impossible to see with a "naked 
eye", they still contribute to the total error budget. For 
this reason we add an additional error in quadrature to 
obtain a reduced equal to 1 for a simple Keplerian 
model. The RVs are flrst modeled with a Keplerian RV 
model and a least-squares orbital flt is made simultane- 
ously to the velocities of the primaries and secondaries 
(the orbital solutions and RVs will be published sepa- 
rately). The residuals are then inspected for planetary 
signatures which would obviously be the same for a pri- 
mary and secondary. In the process, we have not found 
any clear planetary signatures. 

The procedure of iCumming et aTl (|1999[ ) for evaluat- 
ing the regions in Mass-Period space in which compan- 
ions can be ruled out was modified for application to the 
fit residuals to the 2-body Keplerian orbits of the target 
binaries in order to search for and place limits on ad- 
ditional components as follows. This implementation of 
the analysis algorithm w as thoroughly teste d during the 
SIM Double-Blind-Test (|Traub et al.ll2009al lbh. 

A grid of potential companion orbital periods is sam- 
pled logarithmically at values of P = 2fT/I, where / 
is an optional excess factor (here 4) for finer sampling, 
T is the timespan of the observations, and / is an inte- 
ger beginning at 1 and continuing through that at which 
the sampled period is ~ 5x the binary orbital period 
(at which point many companion orbits are dynamically 
unstable). The exact limiting orbital periods shown for 
each star below were cal c ulated using the equation (3) 
from iHolman fc Wiegertl (|1999f ) . At each sampled or- 
bital period, the RV residuals for each star are fit to a 
Keplerian orbit with floating offsets as: 



V = Acos2TTt/P + Bsm2TTt/P - 



- Vo,, 



(5) 



where Vo^i is a floating velocity offset (the i represent- 
ing that different floating offsets are used for each star 
and each observatory; e.g. for binary stars with Keck 
data from before and after the detector upgrade, TNG- 
SARG data, and Lick data, a total of 8 different veloc- 
ity offset parameters are used). At each sample period, 
the x^-niinimizing values of A, B, and Wo,i and result- 
ing value are evaluated. The circular orbit model is 
fo und to be sensitive to m ildly eccentric companions, as 
in ICumming et ahl ()1999l ). 

The values of the best-fit circular orbits were 
recorded for each potential companion period. These 
were converted to the z-statistic of ICumming et al.l 
(|1999f ) for plotting signal power in a periodogram as: 



Q 



Xo 



v^min, circ 



(6) 
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0,„ = ^^£^(^ 

where is twice the number of double-RV measure- 
ments, Ny^ . is the number of independent velocity off- 
sets Vo.i, Xo is the fit to no orbit at all {v = Vo,i)- (When 
making comparisons for calculating false alarm probabil- 
ities and detection limits, the Q coefficients cancel and 
their actual values have little impact.) 

Once the maximum value of z has been found for the 
entire data set, the false alarm probability (FAP) of that 
largest signal is calculated as follows. A Gaussian ran- 
dom number generator is used to create 10,000 synthetic 
data sets with no Keplerian signal but with the same 
cadence and measurement uncertainties as the original 
data. T his introduces a s light difference with the proce- 
dure of iCumming et all ()1999( ): a random number gen- 
erator is used instead of scrambling and rescaling the 
fit residuals. As noted by those authors, the difference 
in these approaches causes little change in the results. 
These synthetic zero-signal data sets are analyzed with 
the same procedure as for the original data residuals. 
Because these synthetic data sets are known to contain 
no real signal, the fraction of these whose analysis show 
power larger than the maximum z of the real data set 
determines how likely that most significant value is to be 
a false positive. 

Finally, the threshold companion RV signal (as a func- 
tion of orbital period) that can be excluded using the 
current data residuals is evaluated. At each sample pe- 
riod, an initial guess for the RV amplitude K that can be 
excluded is made, and 10,000 synthetic data sets with a 
signal of that amplitude are generated. The orbital phase 
is selected randomly with flat distribution. The fraction 
of orbits with z exceeding that of the data is computed, 
the RV amplitude K is modified, and the procedure it- 
erated until the fraction is constrained near the desired 
reliability value (here 99%). In Figures 6-9 (right pan- 
els), the z periodograms are shown for each binary, with a 
horizontal line at the 1% FAP level as determined from 
the signal-free synthetic data sets. In Figures 6-9 (left 
panels), the Mass(sin i)-Period phase-space in which cir- 
cumbinary companions can be excluded at the 99% confi- 
dence level are shown, where the values of Mass(sin i) are 
calculated from the threshold companion RV amplitudes 
K using the total masses of the binary, as estimated for 
each system. We have also tested orbits with non-zero 
eccentricities. It turns out that moderate eccentricities 
make little impact as the threshold lines move up by a 
factor of about 1.5. This is demonstrated on the case of 
HD9939 in Figure 10. 

The targets and their RV data sets are summarized in 
Table 1 where Sp denotes the target's spectral type (ei- 
ther combined or for each component). Mi. 2 the masses 
of the primary and secondary used in the analysis, 
^orbital the Orbital period, ^stable and ^stable sre the 
semi-major axis and orbital period of the first stable or- 
bit of a circumbinary planet computed using the equa- 
tion (3) from Holm an fc Wie gert (1999). Ni^2 denotes 
the number of available RV measurements for the entire 
data set and each subset, rmsi^2 the corresponding best- 
fit rms from a Keplerian orbit for the entire data set as 
well as subsets, Tspan the time span of the entire RV 
data set. ai^2 are the formal RV errors and ei^2 are the 



additional errors added in quadrature. 

5. CONCLUSIONS 

Our novel iodine cell based RV technique allows one to 
measure precise RVs of the components of double-lined 
spectroscopic binary stars. With this technique, in 2003 
we have initiated TATOOINE, a radial velocity search for 
circumbinary planets around a sample of ~50 SB2s. In 
this paper we present the first results from this survey — 
a non detection of exoplanets in the 0.3-3 Mj^p regime 
with the orbital periods of up to 5.3 years around 10 
SB2s. 

Recently two circ umbinary plane ts around an eclips- 
ing binary HW Vir ()Lee et al.ll2009[ ) and a circumbinary 
brown dwar f arou nd en eclipsing binary HS0705-I-6700 
(jOian et al.l l2009[ ) have been claimed to be detected 
by means of eclipse timing. This is however not the 
first time when substellar companions or planets have 
been detected with a timing technique. In addition to 
the confirmed case of the three rocky planets around 
a mi llisecond pulsar B1257-I-12 (|Konacki fc WolszczanI 
'20031), planets have been claimed to orbit a pulsar 
B0329-I-54 based on the tiniing of its radio pulses 
(jPemianski fc Proszynskiiri979l : lShabanoval[l995f l. Later 
it was demonstrated that the timi ng variation is quasi- 
periodic and is not due to planets (|Konacki et aLl flgggi. 
One is left to wonder if the two cases of the eclipse tim- 
ing variations are indeed best explained by a periodic 
signal due to circumbinary bodies and not an unrecog- 
nized quasi-periodic phenomenon mimicking a periodic 
planetary signal. 

Our sample is too small to risk any decisive conclusion 
about the frequency of circumbinary planets. This re- 
minds us of the pioneering search for planets around sin- 
gle stars by Campbell et al. (1988). In particular, their 
non-detection of the population of hot Jupiters. We also 
have not detected any short period planets or to be pre- 
cise planets with periods near the inner orbital stability 
limits. The question of the existence of such circumbi- 
nary planets and hence the impact of the central body's 
binarity on the parking mechanisms for migrating planets 
remains open. It should also be noted that even though 
all our targets are SB2s, they constitute quite a diverse 
sample. Their orbital periods range from 5.7 to 57 days 
and the primary to secondary mass ratios range from 1 
to 1.6. This presumably should make an impact on the 
formation and evolution of circumbinary planets. The 
targets from RV surveys for planets around single stars 
are more homogenous and provide for an overall similar 
environment. 

When comparing our planet detection limi ts to those 
from the Anglo- Australian Planet Search (O' Toole et al.l 
'2009), th e McDonlad Obseravtory Planet Search 

t Wittenmver et al.l l2006f) and the Keck Planet Search 
Gumming et al.l 120081 ) fall employing iodine cells), one 
will notice that these surveys allow for a detection of 
up to several times less massive planets. This is due to 
a higher precision (now approaching 1 ms~^ in many 
cases) , typically larger number of RV measurements and 
a longer time span of the data sets. On the other hand, 
the comparison wit h the planet detection limits for the 
Lick Planet Search (jCumming et al.l 119990 . the first RV 
survey to employ an iodine cell, demonstrates that our 
survey performs similarly. One should also remember 
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that the mass of the central body in our survey is typi- 
cally 1.5-2.0 times higher then in the case of single stars 
which obviously decreases the sensitivity to planets in 
terms of masses by such a factor. 

Our search continues and the changes to the method 
are constantly being made to improve the RV precision. 
Theoretical works on the formation and evolution of cir- 
cumbinary planets are strongly encouraged as now we 
have the observational tools to test them. 

We thank Lucasfilm Ltd for inspiring the TATOOINE 
planet search (and careers of many of us), and Lu- 
casfilm's Senior Director of Business Affairs David An- 
derman for an excellent tour of the Lucasfilm complex 
upon hearing about our program. The tour of Lu- 
casfilm was a highlight of the undergraduate research 
experiences of Agnieszka Czeszumska, Sam Halverson, 
Tony Mercer, and Jackie Schwehr. We thank the Cal- 
ifornia and Carnegie Exoplanet Search team, and Ge- 



off Marcy in particular, for allowing us access to their 
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ated at CDS, Strasbourg, France. MWM acknowledges 
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Double Blind Test, under JPL contract 1336910. This 
research has made use of the Simbad database, oper- 
ated at CDS, Strasbourg, France. The observations on 
the TNG/SARG have been funded by the Optical In- 
frared Coordination network (OPTICON), a major in- 
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frastructures Programme of the European Commissions 
Sixth Framework Programme. 
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Fig. 1. — The RVs of HD9939 (top) and HD13974 (bottom) as a function of the orbital phase (a), and the residuals (observed minus 
modeled RVs) as a function of the orbital phase (b) and time (c). The primary is denoted with filled symbols, the secondary with open ones 
and the best-fit RV model with a solid line. The Keck I/Ifires is denoted with circles, Shane/CAT/Hamspec with triangles and TNG/Sarg 
with stars. 
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Note. — HD13974 is somewhat challenging to disentangle as the RV amplitudes are only 10.2 and 14.2 kms~^ for the primary and secondary respectively. Its mass seems to be low 
by a factor of 2, but we are using this published value anyways. HD221950 has relatively wide spectral lines compared to the remaining targets in this sample and this presumably has 
an impact on the tomographic disentangling and the final RV precision. 



